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We establish the cutoff-dependent upper Higgs boson mass bound by means of direct lattice 
computations in the framework of a chirally invariant lattice Higgs- Yukawa model emulating the 
same chiral Yukawa coupling structure as in the Higgs-fermion sector of the Standard Model. As 
expected from the triviality picture of the Higgs sector, we observe the upper mass bound to decrease 
with rising cutoff parameter A. Moreover, the strength of the fermionic contribution to the upper 
mass bound is explored by comparing to the corresponding analysis in the pure $ 4 -theory. Our final 
results on the cutoff-dependent upper Higgs boson mass bound are summarized in Fig. \T\ 
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FIG. 1: The cutoff dependence of the upper Higgs boson mass bound is presented in panel (a) as obtained 
from the infinite volume extrapolation results in Tab. [V] The dashed and solid curves are fits of the data 
arising from the full Higgs- Yukawa model (HY) and the pure <3? 4 -theory, respectively, with the analytically 
expected cutoff dependence in Eq. (|47)l. Panel (b) shows the aforementioned fit curves extrapolated to larger 
values of the cutoff A. In both panels the highlighted bands reflect the uncertainty of the respective fit curves. 



I. INTRODUCTION 

Given the existing evidence for the triviality of the Higgs sector [l|-0j of the Standard Model, the 
latter theory can only be considered as an effective description of Nature valid at most up to some 
cutoff scale A. The Higgs sector is thus intrinsically connected with a finite, but unknown cutoff 
parameter A that cannot be removed. Beyond that threshold an extension of the theory will finally 
be required. Apriori, the size of this scale A, at which the Standard Model would need such an 
extension, is unspecified. However, the potential discovery of the Higgs boson at the LHC (as well as 
its non-discovery together with corresponding exclusion limits) can shed light on this open question. 
This can, for instance, be achieved by comparing the experimentally revealed Higgs boson mass or its 
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exclusion limits, respectively, with the cutoff-dependent upper and lower Higgs boson mass bounds 
arising in the Higgs sector of the Standard Model. 

Besides the obvious interest in narrowing the interval of possible Higgs boson masses consistent 
with phenomenology, the latter observation was the main motivation for the great efforts spent on the 
determination of cutoff-dependent upper and lower Higgs boson mass bounds. In perturbation theory 
such bounds have been derived from the criterion of the Landau pole being situated beyond the cutoff 
of the theory (see e.g. @Jl3Q, from unitarity requirements (see e.g. [IiT[l3l ]) and from vacuum stability 
considerations (see e.g. [isl. Il4l - [l8j ). as reviewed in Ref. [l9j . 

However, the validity of the perturbatively obtained upper Higgs boson mass bounds is unclear, 
since the corresponding perturbative calculations had to be performed at rather large values of the 
renormalized quartic coupling constant. The latter remark thus makes the upper Higgs boson mass 
bound determination an interesting subject for non-perturbative investigations, such as the lattice 
approach. 

The main objective of lattice studies of the pure Higgs and Higgs- Yukawa sector of the electroweak 
Standard Model has therefore been the non-perturbative determination of the cutoff dependence of 
the upper Higgs boson mass bounds There are two main developments that warrant the 

reconsideration of these questions. First, with the advent of the LHC, we are to expect that the mass 
of the Standard Model Higgs boson, if it exists, will be revealed experimentally. Second, there is, in 
contrast to the situation of earlier investigations of lattice Higgs- Yukawa models [2a - |29j , which suffered 
from their inability to restore chiral symmetry in the continuum limit while lifting the unwanted 
fermion doublers at the same time, a consistent formulation of a Higgs- Yukawa model with an exact 
lattice chiral symmetry [30j based on the Ginsparg- Wilson relation [3l| . This new development allows 
to maintain the chiral character of the Higgs-fermion coupling structure of the Standard Model on 
the lattice while simultaneously lifting the fermion doublers, thus eliminating manifestly the main 
objection to the earlier investigations. The interest in lattice Higgs- Yukawa models has therefore 
recently been renewed (32T - [3l^ | . In particular, the phase diagram of the new, chirally invariant Higgs- 
Yukawa model has been discussed analytically by means of a large Nf calculation [3^, H3| as well as 
numerically by direct Monte-Carlo computations [36| . Moreover, the lower Higgs boson mass bounds 
derived in this lattice model have been presented in Ref. [38[ . A comprehensive review of these results 
can be found in Ref. [39[ . 

In the present paper we intend to determine the dependence of the upper Higgs boson mass bound 
on the cutoff parameter A by direct Monte-Carlo calculations. In sections [H] and IIIII we begin this 
venture by introducing the considered chirally invariant lattice Higgs- Yukawa model and discussing 
the actual simulation strategy, respectively. Details about the determination of the properties of the 
Goldstone and Higgs boson, in particular their renormalized masses, are then given in sections IIVI 
and [V] As a crucial step towards the final determination of the upper mass bound we confirm in 
section [VI] that the largest Higgs boson masses are indeed obtained at infinite bare quartic coupling, 
as expected from perturbation theory. We then present our results on the cutoff dependence of the 
upper Higgs boson mass bound in section fVIII and examine also the encountered finite volume effects. 
Eventually, the lattice data on the Higgs boson mass bounds are extrapolated to the infinite volume 
limit, yielding then our final result already presented in Fig. [TJ 



The model that will be considered in the following, is a four-dimensional, chirally invariant SU(2) l x 
U(1)y lattice Higgs- Yukawa model based on the Neuberger overlap operator [io, Sl|, aiming at the 
implementation of the chiral Higgs-fermion coupling structure of the pure Higgs- Yukawa sector of the 
Standard Model reading 



with (p = iTitp* , Tj being the Pauli matrices, and yt^b denoting the bare top and bottom Yukawa 
coupling constants. In this model the consideration is restricted to the top-bottom doublet (i, b) 
interacting with the complex Higgs doublet ip, which is a reasonable simplification, since the Higgs 
dynamics is dominated by the coupling to the heaviest fermions (apart from its self-coupling). 



II. THE SU(2) i x U(l) y LATTICE HIGGS- YUKAWA MODEL 




(1) 
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The fields contained within the lattice model are thus the scalar field ip, encoded here however in 
terms of the four-component, real scalar field $ for the purpose of a convenient lattice notation, as well 
as Nf top-bottom doublets represented by eight-component spinors tp® = (tw,&W), i — l,...,Nf. In 
this approach the chiral character of the targeted coupling structure in Eq. ([I]) will be preserved on 
the lattice by constructing the fermionic action Sf from the Neuberger overlap operator 2?( 0-u ) acting 
on the aforementioned fermion doublets. The overlap operator is given as 

= Jl + -A=}, A = VW- P , 0<p<2r (2) 

where p is a free, dimensionless parameter within the specified constraints that determines the radius 
of the circle formed by the entirety of all eigenvalues of T>( ov ^ in the complex plane. The operator 
T)( w > denotes here the Wilson Dirac operator defined as 

-D m =Y,^-^>l ( 3 ) 

where V^ :b ' s are the forward, backward and symmetrized lattice nearest neighbor difference operators 
in direction p, while the so-called Wilson parameter r is chosen here to be r = 1 as usual. 

The overlap operator was proven to be local in a field theoretical sense also in the presence of 
QCD gauge fields at least if the latter fields obey certain smoothness conditions [42l.l43j . The locality 
properties were found to depend on the parameter p and the strength of the gauge coupling constant. 
At vanishing gauge coupling the most local operator was shown to be obtained at p — 1. Here, the 
notion 'most local' has to be understood in the sense of the most rapid exponential decrease with 
the distance \x — y\ of the coupling strength induced by the matrix elements r D < £$ between the field 
variables at two remote space-time points x and y. For that reason the setting p = 1 will be adopted 
for the rest of this work. 



Exploiting the Ginsparg- Wilson relation [31| as proposed in Ref. [30| one can then write down a 
chirally invariant SU(2) L x U(l) y lattice Higgs- Yukawa model according to 

Z = { DSD^D^e- 5 *^-^*'^ with (4) 



Nf 

S F [*,il>,4] = Y,$ (i) I? (o, ' ) +P + t diag(y t ,y b )P + +P_diag(y t ,y fc )0P_ (5) 



i=l 



M 



where the particular form of the 0(4)-symmetric purely bosonic action S$[$] will be given later. 
It is further remarked that the four-component scalar field <5> Xl defined at the Euclidean site indices 
x = (t, x) of a x L t -lattice, has been rewritten here as a quaternionic, 2x2 matrix <fr x = ^0^, 9o = 
1, 9j = —iTj with t denoting the vector of Pauli matrices, acting on the flavour index of the fermion 
doublets. The so far unspecified left- and right-handed projection operators P± and their lattice 
modified counterparts P± associated to the Neuberger Dirac operator are given as 

p±=^^, p±=^, j 5 = 7 ji--v^y (6) 

The action in Eq. <j5j) now obeys an exact global SU(2)x, x U(l)y lattice chiral symmetry. For VLl e 
SU(2) and e € 1R the action is invariant under the transformation 

iP^u R P+ij + u L n L P^, $ -> ijp + n{ul + QP-U^, (7) 
<l> -> u^l, ->n L tful (8) 

with Ul,r,c/) = exp(ieY) denoting the respective representations of the global U{l)y symmetry group. 
Employing the explicit form of the hypercharge Y being related to the isospin component 23 and the 
electric charge Q according to Y = Q — I3, the above U(1)y matrices can explicitly be parametrized 
as 

/ +ie/6 \ / 2»£/3 \ / p +ie/2 \ 

U L = e+ ie/ 6 )> U *={ e -ie/3 K ^ = ( e -U/2 > (9) 



4 



for the case of the considered top-bottom doublet. For clarification it is remarked that the right- 
handed fields are isospin singlets and have only been written here in form of doublets for the sake of a 
shorter notation. Note also that in the mass-degenerate case, i.e. y t = yb, the above global symmetry 
is extended to SU(2)l x SU(2)_r. In the continuum limit the modified projectors P± converge to P± 
and the symmetry in Eq. (J7HSJ) thus recovers the continuum SU(2)i x U(l)y global chiral symmetry 
such that the lattice Higgs- Yukawa coupling becomes equivalent to Eq. (JTJ when identifying 

V* = C ■ cs§_ ofS > ?* = WP* = C -\ _*2 , s li , Vt,b = —pr (10) 



for some real, non-zero constant C . 

The so far unspecified purely bosonic action S$> is chosen here to be the lattice version of the $ 4 - 
action parametrized in terms of the hopping parameter n and the lattice quartic coupling constant A 
according to 

5$ = - K e *t + + E $ ^ + A E - N f) 2 ' ( n ) 

X X 

which is a convenient parametrization for the actual numerical computations. However, this form of 
the lattice action is fully equivalent to the lattice action in continuum notation 

S M = E^ v ^- v ^-+E^ m o^^ + E A K^) 2 ' ( 12 ) 

X.f^i X X 

given in terms of the bare mass mo, the bare quartic coupling constant A, and the lattice derivative 
operator Vjj. The aforementioned connection can be established through a rescaling of the scalar field 
$ and the involved coupling constants according to 

rr-f^l+i^W , A 2 1 - 2N } \ - 8k y t , b 

Finally, the potential appearance of a sign problem in the framework of the introduced Higgs- 
Yukawa model shall be briefly addressed. In the mass-degenerate case, i.e. for y t = yb, one finds that 
det(A4) £ K, since all eigenvalues of Ai come in complex conjugate pairs according to 

VMV*=M*, with V = 7o7275T 2 . (14) 

This is in contrast to the general case with yt ^ yb, where the above relation no longer holds. Through- 
out this work we will therefore only consider the aforementioned mass-degenerate scenario, where the 
top and bottom quarks are assumed to have equal masses, to certainly exclude any complex-valued 
phase of the fermion determinant. This, however, still leaves open the possibility of an alternating 
sign of det(Ai). We have therefore explicitly monitored the sign of det(Ai) but did never encounter 
any sign alteration in our actually performed Monte-Carlo computations, meaning that the numerical 
calculations in the mass-degenerate case are perfectly sane. A more detailed discussion of the phase 
of the fermion determinant in the non-degenerate case can be found in Ref. [39| . 



III. SIMULATION STRATEGY AND CONSIDERED OBSERVABLES 



The eventual aim of this work is the non-perturbative determination of the cutoff-dependent upper 
bound of the Higgs boson mass. The general strategy that will be applied for that purpose is to 
scan through the whole space of bare model parameters searching for the largest Higgs boson mass 
attainable within the pure Higgs- Yukawa sector at a fixed value of the cutoff, while being in consistency 
with phenomenology. This will be done by numerically evaluating the finite lattice model of the Higgs- 
Yukawa sector introduced in the preceding section and extrapolating the obtained results to the infinite 
volume limit. 

The crucial idea is that the aforementioned requirement of reproducing phenomenology restricts 
the freedom in the choice of the bare model parameters itIq, yt,b, A. For that purpose we exploit 
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here the phenomenological knowledge of the renormalized quark masses and the renormalized vacuum 
expectation value of the scalar field (vev). For the reasons given in the previous section, however, 
the top and bottom quarks will be considered to be mass-degenerate. Throughout this work m t /a = 
rah I a — 175 GeV and v r /a — 246 GeV will be assumed. Here m*, mj, and v r are the renormalized 
top and bottom quark masses as well as the renormalized vev in dimensionless lattice units, while 
a denotes the lattice spacing. The aforementioned three conditions leave open an one-dimensional 
freedom in the bare parameters, which can be parametrized in terms of the bare quartic self-coupling 
constant A. However, aiming at the upper Higgs boson mass bounds, this remaining freedom can be 
fixed, since it is expected from perturbation theory that the lightest Higgs boson masses are obtained 
at vanishing self-coupling constant A = 0, while the heaviest masses are attained at infinite coupling 
constant A = oo, respectively. That this conjecture actually holds also in the non-perturbative regime 
of the model, i.e. at large values of A, is explicitly demonstrated in section fVTl allowing then to restrict 
the search for the upper mass bound to the setting A = oo. 

Furthermore, the model has to be evaluated in the broken phase, i.e. at (ip) ^ 0, to respect the 
observation of spontaneous symmetry breaking, however close to a second order phase transition to 
a symmetric phase to allow for arbitrarily large correlation lengths as required in any attempt of 
pushing the cutoff parameter to arbitrarily large values. 

However, in the given lattice model the expectation value (ip) would always be identical to zero 
due to the symmetries in Eq. J7][H]). The problem is that the lattice averages over all ground states 
of the theory, not only over that one which Nature has selected in the broken phase. To study the 
mechanism of spontaneous symmetry breaking nevertheless, one usually introduces an external current 
J, selecting then only one particular ground state. This current is finally removed after taking the 
thermodynamic limit, leading then to the existence of symmetric and broken phases with respect to 
the order parameter (ip) as desired. An alternative approach, which was shown to be equivalent in the 
thermodynamic limit [44j-|46j, is to define the vacuum expectation value (vev) v as the expectation 
value of the rotated field ip rot given by a global transformation of the original field <p according to 

p:f = u[p\p x (is) 

with the SU(2) matrix U[<p] selected for each configuration of field variables {tp x } such that 





E 



,rot 
'fix 



(16) 



Here we use this second approach. According to the notation in Eq. (|12p . which already includes a 
factor 1/2, the relation between the vev v and the expectation value of ip rot is then given as 

(O = ( I ) • (17) 

In this setup the unrenormalized Higgs mode h x and the Goldstone modes gl,g x , g x , can then directly 
be read out of the rotated field according to 

..rot _ ( 9x + ^9x | / 1 o\ 

-\v + h x -igl)- ^ 

The great advantage of this approach is that no limit procedure J —¥ has to be performed, which 
simplifies the numerical evaluation of the model tremendously. 

The physical scale of the lattice computation, i.e. the inverse lattice spacing a -1 , can then be 
determined by comparing the renormalized vev v r = v j \J Zq measured on the lattice with its phe- 
nomcnologically known value according to 

246 GeV = - = _JL , (19) 



where Zq denotes the Goldstone renormalization constant. The cutoff parameter A of the underlying 
lattice regularization, which is directly associated to the lattice spacing a, can then be defined as 



A = a -\ (20) 



G 



Of course, this definition is not unique and other authors use different definitions, for instance 
A = ir/a motivated by the value of the momenta at the edge of the Brillouin zone. However, since 
the quantities that actually enter any lattice calculation are rather the lattice momenta p^ — sm(p M ) 
instead of the momenta p^,, which are connected through the application of a sine function, it seems 
natural to choose the definition of the cutoff A given in Eq. (j2T)f . 

Next, the extraction technique for the Goldstone renormalization constant entering Eq. (|19[) needs 
to be determined. In the Euclidean continuum the Goldstone and Higgs renormalization constants, 
more precisely their inverse values Zq 1 and Z^ 1 , are usually defined as the real part of the derivative 
of the inverse Goldstone and Higgs propagators in momentum space with respect to the continuous 
squared momentum p 2 at some scale p 2 — — /ig and p 2 = —fJ,jj, respectively. The restriction to the 
real part is introduced to make this definition applicable also in the case of an unstable Higgs boson, 
where the massless Goldstone modes induce a branch cut with discontinuous complex contributions 
to the propagator at negative values of p 2 . This is the targeted definition that shall also be adopted 
to the later lattice calculations. 

On the lattice, however, the propagators are only defined at the discrete lattice momenta p^ = 
2Tm^/L s j, n M = 0, . . . ,L Sjt - 1 according to 



Gh(j>) = (hph-p), 
1 3 



(21) 
(22) 



Q = l 



where the Higgs and Goldstone fields in momentum representation read 



(23) 



with V — L z s ■ Lt denoting the lattice volume. 

Computing the derivative of the lattice propagators is thus not a well-defined operation. Moreover, 
the lattice propagators are not even functions of p 2 , since rotational invariance is explicitly broken 
by the discrete lattice structure. To adopt the above described concept to the lattice nevertheless, 
some lattice scheme has to be introduced that converges to the continuum definitions of Zq and Z^ 
in the limit a — » 0. Here, the idea is to use some analytical fit formulas fcip), Jh{p) derived from 
rcnormalizcd perturbation theory in the Euclidean continuum to approximate the measured lattice 
propagators Gc{p) and Gh{p) at small momenta p 2 < 7 (with p 2 ^ = 4sin 2 (p p /2)) for some appropriate 
value of 7 such that the discretization errors are acceptable. The details of this fit procedure are 
discussed in sections IIVI and |V] One can then define the analytically continued lattice propagators as 



G c G (Pc) = fc(Pc) and G c H {p c ) = f H (p c )- 



(24) 



In the on-shell scheme the targeted Goldstone and Higgs renormalization constants Zq and Zh can 
then be defined (implicitly assuming an appropriate mapping p c ■<-> p 2 ) as 



to = ( G G ( P i 



_d_ 

:lp°- 
d 
dp 2 



■Re 



g c h {pI) 



(25) 
(26) 



with /ig. = rriQ and ^i 2 H = m 2 H , where the underlying physical masses mc, mjj are given by the poles 
of the respective propagators on the second Riemann sheet. To adopt this definition to the introduced 
lattice scheme we define the Higgs boson mass run, its decay width Th, and the mass mc of the 
stable Goldstone bosons through 



G C H 11 {im H +r H /2, 0,0,0) =0, and G G (im G , 0, 0, 0) =0 



(27) 



where G C H jj(p c ) denotes the analytical continuation of G c H (p c ) onto the second Riemann sheet. 
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Extracting the Higgs boson mass ran and its decay width Th from simulation data according to 
this definition would, however, require an explicit analytical continuation of the Higgs propagator onto 
the second Riemann sheet, which is beyond our ambitions in this study. 

Following the proposal in Ref. Q the Goldstone and Higgs renormalization factors are rather de- 



termined at the scales fx G 



and fj,jj = m Hp given by the masses and mc P , which will be 



referred to in the following as propagator masses in contrast to the pole masses m# and mo- We thus 
define 

Z G = Z G (m 2 Gp ) and Z H = Z H {m 2 Hp ) 1 (28) 

where the propagator masses m_f/p, mcp are defined through a vanishing real-part of the inverse 
propagators according to 



Re 



and Re 



Gp 



0. 



Pf: ~ 



(29) 



The reasoning for selecting these latter definitions of the Higgs and Goldstone masses is, that the 
required analytical continuation in the case of the Higgs propagator is much more robust, since it only 
needs to extend the measured lattice propagator to purely negative values of p 2 in contrast to the 
situation resulting from the definition in Eq. (|27[) . It is remarked here that the Goldstone propagator 
mass map was only introduced for the sake of an uniform notation, since vhq is identical to map, due 
to the Goldstone bosons being stable particles. 

As for the unstable Higgs boson, however, one finds that the discrepancy between the pole mass 
vriH and the propagator mass mnp is directly related to the size of the decay width Th- In the 
weak coupling regime of the theory the two mass definitions mn and mup can thus be considered to 
coincide up to small perturbative corrections, due to a vanishing decay width in that limit. For the 
pure $ 4 -theory the deviation between mnp and tuh has explicitly been worked out in renormalized 
perturbation theory Q . In infinite volume the finding is 



m H = m Hp ■ I 1 + — (n - l) 2 



4! ■ A r 
16tt 2 



-i 2 



•0(A?) 



(30) 



where A r denotes the renormalized quartic self-coupling constant and n is the number of components 
of the scalar field $, i.e. n — A for the here considered case. This calculation was performed in the 
pure <I> 4 -theory, thus neglecting any fermionic degrees of freedom, and for exactly massless Goldstone 
particles. However, one learns from this result that the definition of mnp in Eq. (I29p as the Higgs boson 
mass is very reasonable at least for sufficiently small values of the renormalized coupling constants. 
The actual discrepancy between ran and mnp as obtained by direct lattice computations of their 
respective definitions in Eq. (|27|) and Eq. (|29|) will explicitly be examined in section [V] for some 
physically relevant parameter setups. It will then indeed be found to be negligible with respect to the 
reachable statistical accuracy. 

The definition of the renormalized quartic self-coupling constant X r that was used in the derivation 
of Eq. d30j) is 



rr> 



Hp 



m 



Gp 



8v 2 



(31) 



which shall also be taken over to the considered Higgs- Yukawa model. In principle, it would also 
be possible to determine the renormalized quartic coupling constant X r through the evaluation of 
the amputated, connected, one-particle-irreducible four-point function at a specified momentum con- 
figuration as it is usually done in perturbation theory. However, the signal to noise ratio of the 
corresponding lattice observable is suppressed by the lattice volume. It is thus extremely hard to 
measure the renormalized quartic cou plin g constant in lattice calculations by means of the direct 
evaluation of such four-point functions [47] . Instead, the alternative definition of A r given in Eq. pip 
will be adopted here. It is further remarked that this definition was shown Q to coincide with the 
bare coupling parameter A to lowest order in the pure 4> -theory. 

Regarding the top and bottom quark fields, we are here only interested in the corresponding masses 
mt,TOfe. These can directly be obtained by studying the fermionic time correlation functions C/(At) 
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at large Euclidean time separations At, where f — t,b denotes the quark flavour here. On the lattice 
the latter time correlation functions can be defined as 

C f (At) = j—js E ( 2ReTr (/£.*+At,a ' fR,t,y) ), (32) 

4 s t=0 x,y 

where the left- and right-handed spinors are given through the projection operators according to 

=^-(ft) and (*,&)* = (*,5)P_. (33) 

It is remarked that the given fermionic correlation function would be identical to zero due to the exact 
lattice chiral symmetry obeyed by the considered Higgs- Yukawa model, if one would not rotate the 
scalar field <p according to Eq. (JTSJ) as discussed above. This rotation is implicitly assumed in the 
following. Furthermore, it is pointed out that the full all-to-all correlator as defined in Eq. (|32l) can 
be trivially computed. This all-to-all correlator yields very clean signals for the top and bottom quark 
mass determination. 

The lacking definition of the renormalized Yukawa coupling constants can now be provided as 

ffH m b 

Vt.r = — and y b . r = — , (34) 

v r v r 

reproducing the bare Yukawa coupling constants ytj> at lowest order. According to the presented 
simulation strategy the aim would thus be to tune the above renormalized Yukawa coupling constants 
such that their physically known values would be reproduced in the actual lattice computations. 
However, for having some initial guess for the latter adjustment at hand the tree-level relation 

Vt.b = (35) 

v r 

will be used throughout this work to set the bare Yukawa coupling constants in the lattice computa- 
tions. Comparing the physical fermion masses actually generated in these lattice calculations with the 
targeted ones would then allow to fine tune the Yukawa coupling constants in an iterative refinement 
approach. However, it turns out that this tree-level fixation ansatz already yields quite satisfactory 
results regarding the discrepancy between the targeted and the actually observed quark masses with 
respect to the reached statistical accuracy. 



IV. ANALYSIS OF THE GOLDSTONE PROPAGATOR 



The Goldstone renormalization constant Zq is required for determining the renormalized vacuum 
expectation value v r of the scalar field. It is thus needed for the fixation of the physical scale a -1 
of a given Monte-Carlo run according to Eq. (|19p . This renormalization constant has been defined 
in Eq. (f28|) through a derivative of the inverse Goldstone propagator. As already pointed out in 
section IHlI computing this derivative requires an analytical continuation G c G (p c ) of the discrete lattice 
propagator, which was proposed to be obtained via a fit of the discrete lattice data. 



9 \ 



h \ 



/ h \ 



FIG. 2: Illustration of the diagrams that contribute to the continuous space-time Goldstone propagator Gq(Pc) 
in the Euclidean pure $ 4 -theory at one-loop order. 
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pl \ Vq J 2 Pc \ m G 

4 ■ arctanh ( vl + ~ m ° | with 



The idea here is to construct an appropriate fit function f G (p) based on a perturbative calculation 
of the Goldstone propagator G G (p c ) in continuous Euclidean space-time. In this study the aforemen- 
tioned fit function will only play the role of an effective description of the numerical data to allow for 
the necessary analytical continuation. For its construction we can therefore impose a set of simplifi- 
cations. In particular, we restrict the consideration here to the pure $ 4 -theory. The reasoning behind 
this simplification is that the purely bosonic four-point interaction is expected to yield the dominant 
contributions to the Goldstone propagator in the targeted strong coupling regime with infinite bare A 
but only moderate values of the bare Yukawa coupling constants. To one-loop order the only momen- 
tum dependent contribution to the Goldstone propagator is thus given by the mixed Higgs- Goldstone 
loop illustrated on the right-hand side of Fig. [3J where the system has been assumed to be in the 
broken phase, as desired. At one-loop order the result for the renormalized Goldstone propagator 
then reads 

G G \p c ) = pl + m% + 8n- 2 \lvZ-[X(pl,m%,m%)-X(-m%,m%,m 2 G )] (36) 
where the one- loop contribution X(p 2 , tti 2 Hi m G ) is given as 

2X( P l m%, m%) = ^. arctanh ( Pj±^_^k) + !^_^ . log ( (3?) 

a 

Pl V VQ 

q = (m G -m 2 H +p 2 c ) 2 +im 2 H p 2 c . (38) 

Concerning the singularities of this expression it is noteworthy to add that the given formula can be 
shown to be finite at p c = for m G =/= 0, as desired. 

In principle, one can directly employ the expression in Eq. (|36p as the sought-after fit function 
f G 1 (p). For clarification it is remarked at this point that instead of fitting the lattice propagator G G (p) 
itself with f G (p), it is always the inverse propagator G^ 1 ^) that is fitted with /^(p) = 1//g(p) i n 
the following. However, for the actual fit procedure of the lattice data a modified version of Eq. (|36l) 
is used given as 

/-i(pS) = P 2 + fh l+ A - [X(p 2 ,m 2 H , m 2 G ) - 1(0, fh 2 H , fh G )} 

where an appropriate mapping p 2 ■<-> p is implicit and A, Zq, fh G , rhn are the free fit parame- 
ters. Two modifications have been applied here to the original result. Firstly, the constant term 
I(—fh G ,fhff,m G ) in Eq. (|36[) has been replaced by 1(0, fh 2 Hl m G ) simply for convenience. Since the 
Goldstone mass is close to zero anyhow, this simplification is insignificant for a practical fit procedure. 
For clarification it is recalled that in the presented approach the Goldstone mass m G is actually not 
determined through the nominal value of the latter fit parameter fh G itself, but through the pole 
of the resulting analytical continuation G G (p c ) according to Eq. (|!2T)) . This is also indicated by the 
chosen notation introducing the symbol fh G in addition to the actual Goldstone mass m G . 

More interestingly, however, a global factor Zq has been included in the denominator of Eq. (|39[) in 
the spirit of a renormalization constant. This modification is purely heuristic and its sole purpose is 
to provide an effective description of the so-far neglected fermionic contributions, which is all we need 
at this point. 

Of course, it would be more appropriate to construct a fit ansatz from the renormalized result 
of the Goldstone propagator derived in the full Higgs- Yukawa sector. This would indeed place the 
fit procedure on an even better conceptual footing. However, it will turn out, that the given ansatz 
already works satisfactorily well for our purpose, which is not too surprising, due to the aforementioned 
dominance of the quartic coupling term in that model parameter space being of physical interest here. 

More important seems to be the question what part of the lattice Goldstone propagator G G (p) 
one should actually include into the fit procedure. It was already pointed out in section Mil that 
the consideration of the lattice propagator has to be restricted to small lattice momenta in order 
to suppress contaminations arising from discretization effects. For that purpose the constant 7 has 
been introduced specifying the set of momenta underlying the fit approach according to p 2 < 7. In 
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principle, one would want to choose 7 as small as possible. In practice, however, the fit procedure 
becomes increasingly unstable when lowering the value of 7, because less and less data are then 
included within the fit. In the following example lattice computations, demonstrating the evaluation 
approach for Zq and ma, we will consider the settings 7 = 1, 7 = 2, and 7 = 4. To make the 
discretization effects associated to these not particularly small values of 7 less prominent in the 
intended fit procedure, the inverse lattice propagator G^ijp) is actually fitted with f^ij) 2 ) instead 
of (p 2 ), being a function of the squared lattice momentum p 2 , which is completely justified in the 
limit 7 — >• 0. 



L 3 S x L t 


N f 




A 


Vt 


Vb/yt 


V 


A 


32 4 


1 


0.30039 


00 


0.55139 


1 


0.1008(3) 


2373.0 ± 6.4 GeV 


32 4 


1 


0.30400 


00 


0.55038 


1 


0.1547(1) 


1548.1 ± 1.8 GeV 



TABLE I: The model parameters of the Monte-Carlo runs constituting the testbed for the subsequently 
discussed computation schemes are presented together with the obtained values of the vacuum expectation 
value v and the cutoff A determined by Eq. (|20|) . The degenerate Yukawa coupling constants have been chosen 
here according to the tree-level relation in Eq. 1)351) aiming at the reproduction of the phenomenologically 
known top quark mass. 

The Goldstone propagators obtained in the lattice calculations specified in Tab. Q] are presented in 
Fig. [3] These numerical data of the inverse Goldstone propagator G G l {p) have been fitted with the fit 
formula f^ip 2 ) given in Eq. (jSHJ). One can observe already from the graphical presentation in Fig. [3] 
that the considered fit ansatz /g(p 2 ) describes the numerical data significantly better than the simple 
linear fit formula 

l- G \f) = (40) 

which is additionally considered here for the only purpose of demonstrating the quality of the applied 
fit ansatz fciP 2 )- 

To find an optimal setting for the threshold value 7, the dependence of the fit results on the latter 
parameter is listed in Tab.|TTl where the presented Goldstone mass and the renormalization factor 
Zq have been obtained according to Eq. ((27)) and Eq. (|28l) from the analytical continuation of the 
lattice Goldstone propagator given by G G (p c ) — fa(Pc) and G G (p c ) — Ig{p c ), respectively. 

At first glance one notices that the linear ansatz Ig(p 2 ) yields more stable results than fc(p 2 )- 
These results are, however, inconsistent with themselves when varying the parameter 7. One can 
also observe in Tab. [TT] that the associated average squared residual per degree of freedom y^/dof 
significantly differs from one at the selected values of 7, making apparent that the simple linear fit 
ansatz is not suited for the reliable determination of the Goldstone propagator properties. 





fit ansatz Jg(P 2 ) 


linear fit ansatz 


a(f) 


K 


7 


Zg 


mo 


X 2 /dof 


Z G 


m G 


X 2 /dof 


0.30039 


1.0 


0.9380(107) 


0.027(15) 


1.00 


0.9422(5) 


0.067(2) 


2.61 


0.30039 


2.0 


0.9431(52) 


0.028(11) 


0.81 


0.9507(3) 


0.089(2) 


4.79 


0.30039 


4.0 


0.9457(27) 


0.033(8) 


0.94 


0.9585(2) 


0.114(2) 


6.19 


0.30400 


1.0 


0.9400(90) 


0.029(10) 


1.41 


0.9403(4) 


0.066(1) 


4.40 


0.30400 


2.0 


0.9426(36) 


0.032(7) 


1.07 


0.9476(2) 


0.084(1) 


6.53 


0.30400 


4.0 


0.9478(18) 


0.038(4) 


1.06 


0.9559(1) 


0.111(1) 


9.67 



TABLE II: The results on the Goldstone renormalization factor Zg and the Goldstone mass tug, obtained 
from the fit approaches fc{p 2 ) and Ig{p 2 ) as defined in Eq. (|39[l and Eq. (|40f) . are listed for several settings 
of the parameter 7 together with the corresponding average squared residual per degree of freedom \ 2 /dof 
associated to the respective fit. The underlying Goldstone lattice propagators have been calculated in the 
Monte-Carlo runs specified in Tab. U 
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FIG. 3: The inverse lattice Goldstone propagators calculated in the Monte-Carlo runs specified in Tab. [J 
are presented versus the squared lattice momenta p 2 together with the respective fits obtained from the fit 
approaches /^(p 2 ) in Eq. (|39|l (red solid line) and ^(p 2 ) in Eq. (|40jl (blue dashed line) with 7 = 4.0. From 
left to right the three panel columns display the same data zooming in, however, on the vicinity of the origin 



at p 2 = 0. 



In contrast to that the more elaborate fit ansatz fa(p 2 ) yields much better values of x 2 /dof being 
close to the expected value of one as can be seen in Tab. HU Moreover, the results on the renormal- 
ization constant Zq and the Goldstone mass vhq obtained from this ansatz remain consistent with 
respect to the specified errors when varying the constant 7. In the following the aforementioned 
quantities Zq and ma will therefore always be determined by means of the here presented method 
based on the fit ansatz /g(p 2 ) with a threshold value of 7 = 4, since this setting yields the most stable 
results, while still being consistent with the findings obtained at smaller values of 7. 



V. ANALYSIS OF THE HIGGS PROPAGATOR 



Concerning the analysis of the Higgs propagator we will follow the same strategy as in the previous 
section. Examples of the lattice Higgs propagator as obtained in the Monte-Carlo runs specified in 
Tab. [I] are presented in Fig. 21 These numerical data have been fitted with the ansatz 

t -u«, f+ml + A- [36 (j(p 2 , fh 2 H ) - D H0 ) + 12 (j(f, m|) - D G0 )] 

f H VP) = — 1 7 " " — ) ( 41 ) 

derived from renormalized perturbation theory in the Euclidean pure <I> 4 -theory at one-loop order. 
The restriction to the pure $ 4 -theory is again motivated by the same arguments already discussed in 
the preceding section. In the above formula the one-loop contribution J{p 2 , m 2 ) is defined as 



9 9 . arctanhfo) / p 2 , „ . 

J(p 2 ,m 2 ) = ^, q=J-Ji— -, (42) 

q y 4to z + p z 
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the constants Dhq, Dqo are given as 

D H0 = J(0,m 2 H ) = l, (43) 
D G0 = J(0,m 2 G ) = l, (44) 

and fftjj, A, Zq are the free fit parameters. The parameter fh G is not treated as a free parameter 
here. Instead it is fixed to the value of mc resulting from the analysis of the Goldstone propagator by 
the method described in the previous section. The sole purpose of this approach is to achieve higher 
stability in the considered fit procedure, which otherwise would yield here only unsatisfactory results 
with respect to the associated statistical uncertainties. 

Again, one can observe, however less clearly as compared to the previously discussed examples of 
the case of the Goldstone propagator, that the more elaborate fit ansatz fnip 2 ) describes the lattice 
data more accurately than the simple linear fit approach 

l H \f) = P -^- (45) 

This is better observable in the lower row of Fig.2]than in the upper row, where the differences tend to 
be rather negligible. The reason why the observed differences between the two fit approaches are less 
pronounced here, as compared to the situation in the preceding section, simply is, that the threshold 
value 7 was chosen here to be 7 = 1 which will be motivated below. This setting of 7 is much smaller 
than the value underlying the previously discussed examples of the Goldstone propagators and causes 
the linear fit to come closer to the more elaborate ansatz /.fj(p 2 ). 
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FIG. 4: The inverse lattice Higgs propagators calculated in the Monte-Carlo runs specified in Tab. U are pre- 
sented versus the squared lattice momenta p 2 together with the respective fits obtained from the fit approaches 



fn 



in Eq. (|4T) (red solid line) and l^if) 2 ) in Eq. (|4"5j) (blue dashed line) with 7 = 1.0. From left to right 



the three panel columns display the same data zooming in, however, on the vicinity of the origin at p =0. 



The Higgs propagator mass m# p defined in Eq. (|29|) and the Higgs pole mass m# together with its 
associated decay width Th given by the pole of the propagator on the second Riemann sheet according 
to Eq. ([2"T]) can then be obtained by defining the analytical continuation of the lattice propagator as 
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G c h(Pc) = Ih{Pc) and G c H (p c ) — Ih(p c ), respectively. The results arising from the considered fit 
procedures are listed in Tab. IIIII for several values of the threshold value 7. However, since the linear 
function Ih{Pc) can not exhibit a branch cut structure, the pole mass equals the propagator mass and 
the decay width is identical to zero when applying the linear fit approach. That is the reason why 
only the Higgs boson mass mu is presented in the latter scenario. We further remark that the values 
of Th arising along with the determination of tuh are - as expected - rather unstable due to the 
required analytical continuation of the propagator onto the second Riemann sheet. We therefore do 
not present these numbers here. 

One observes in Tab. IIIII that the Higgs boson masses obtained from the linear fit ansatz Ih(p 2 ) are 
again inconsistent with the respective results obtained at varying values of the threshold parameter 
7, thus rendering this latter approach unsuitable for the description of the Higgs propagator. This 
becomes also manifest through the presented values of the average squared residual per degree of 
freedom \ 2 /dof associated to the linear ansatz, which are clearly off the expected value of one (with 
the exception of the case of 7 = 0.5). 





Fit ansatz /ff(p 2 ) 


Fit ansatz Ih(p 2 ) 


K 


7 


m,H P 


run 


X 2 /dof 




X 2 /dof 


0.30039 


0.5 


0.253(2) 


0.296(83) 


1.17 


0.253(2) 


1.13 


0.30039 


1.0 


0.252(2) 


0.253(2) 


1.20 


0.261(2) 


1.62 


0.30039 


2.0 


0.246(2) 


0.249(2) 


1.09 


0.273(1) 


2.58 


0.30400 


0.5 


0.405(3) 


0.406(3) 


1.43 


0.399(2) 


1.75 


0.30400 


1.0 


0.409(1) 


0.410(1) 


1.16 


0.409(1) 


2.23 


0.30400 


2.0 


0.409(1) 


0.412(1) 


1.27 


0.423(1) 


4.63 



TABLE III: The results on the Higgs propagator mass rriHp and the Higgs pole mass ran obtained from 
the fit approaches fii(p 2 ) in Eq. (|4lT) and Ih(p 2 ) in Eq. (l45l) are listed for several settings of the parameter 
7 together with the corresponding average squared residual per degree of freedom x 2 /dof associated to the 
respective fit. For the linear fit ansatz ln(p 2 ) only the pole mass is presented, since one finds mu v = mn when 
constructing the analytical continuation G c H {pc) through Ih(Pc)- The underlying Higgs lattice propagators 
have been calculated in the Monte-Carlo runs specified in Tab. [I] 

Again the situation is very different in case of the more elaborate fit ansatz fnip 2 ) yielding 
significantly smaller values of x 2 /dof. The presented results on the propagator mass as well 

as the pole mass mn are also in much better agreement with the corresponding values obtained at 
varied threshold parameter 7. Moreover, the values of mn v and thh are consistent with each other 
with respect to the given errors, finally justifying the identification of the Higgs boson mass with the 
propagator mass rriHp- 

From the findings presented in Tab. IIIII one can conclude that selecting the threshold value to be 
7 = 1 for the analysis of the Higgs propagator is a very reasonable choice, which leads to consistent 
and satisfactory results. This is the setting that will be used for the subsequent investigation of the 
upper Higgs boson mass bounds to determine the propagator mass mup- It is further remarked that 
the here chosen value of 7 is much smaller than the value 7 = 4 selected in the preceding section for 
the analysis of the Goldstone propagator. While this large setting worked well in the latter scenario, 
it leads to less consistent results in the here considered case and has therefore been excluded from the 
given presentation. 

VI. DEPENDENCE OF THE HIGGS BOSON MASS ON THE BARE COUPLING 

CONSTANT A 

We now turn to the question whether the largest Higgs boson mass is indeed obtained at infinite 
bare quartic coupling constant for a given set of quark masses and a given cutoff A as one would 
expect from perturbation theory. Since perturbation theory may not be trustworthy in the regime 
of large bare coupling constants, the actual dependence of the Higgs boson mass on the bare quartic 
coupling constant A in the scenario of strong interactions is explicitly checked here by means of direct 
Monte-Carlo calculations. The final answer of what bare coupling constant produces the largest Higgs 
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boson mass will then be taken as input for the investigation of the upper mass bound in the subsequent 
section. 

For this purpose some numerical results on the Higgs propagator mass niH P are plotted versus the 
bare quartic coupling constant A in Fig. The presented data have been obtained for a cutoff that 
was intended to be kept constant at approximately A w 1540 GeV by an appropriate tuning of the 
hopping parameter, while the degenerate Yukawa coupling constants were fixed according to the tree- 
level relation in Eq. (|35j) aiming at the reproduction of the top quark mass. One clearly observes that 
the Higgs boson mass monotonically rises with increasing values of the bare coupling constant A until 
it finally converges to the A = oo result, which is depicted by the horizontal line in the presented plot. 
From this result one can conclude that the largest Higgs boson mass is indeed obtained at infinite 
bare quartic coupling constant, as expected. The forthcoming search for the upper mass bound will 
therefore be restricted to the scenario of A = oo. 
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FIG. 5: The Higgs boson mass mn v and the renormalized quartic coupling constant A r are shown versus 
the bare coupling constant A in panels (a) and (b), respectively. These results have been obtained in direct 
Monte-Carlo calculations on a 16 3 x 32-lattice with the degenerate Yukawa coupling constants fixed according 
to the tree-level relation in Eq. (|35|l aiming at the reproduction of the top quark mass. The hopping parameter 
was tuned with the intention to hold the cutoff constant, while the actually obtained values of A fluctuate 
here between 1504 GeV and 1549 GeV. The horizontal lines depict the corresponding results at infinite bare 
coupling constant A = oo, and the highlighted bands mark the associated statistical uncertainties. 



Furthermore, the corresponding behaviour of the renormalized quartic coupling constant A r as 
defined in Eq. (|31[) is presented in Fig.(5}D. As expected one observes a monotonically rising dependence 
of A r on the bare coupling constant A, eventually converging to the A = oo result depicted by the 
horizontal line. 



VII. RESULTS ON THE UPPER HIGGS BOSON MASS BOUND 



We now turn to the actually intended non-perturbativc determination of the cutoff-dependent upper 
Higgs boson mass bound m^f(A). Given the knowledge about the dependence of the Higgs boson mass 
on the bare quartic self-coupling constant A the search for the desired upper mass bound can safely be 
restricted to the scenario of an infinite bare quartic coupling constant, i.e. X — oo. Moreover, we will 
restrict the investigation here to the mass degenerate case with y t = yb, since the fermion determinant 
det(A4) can be proven to be real in this scenario as discussed in section UT1 

Concerning the cutoff parameters A that are reachable with the intended lattice calculations, a 
couple of restrictions limit the range of the accessible energy scales. On the one hand all particle 
masses have to be small compared to A to avoid unacceptably large cutoff effects, on the other 
hand all masses have to be large compared to the inverse lattice side lengths to bring the finite 
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volume effects to a tolerable level. As a minimal requirement we demand here that all particle masses 
rh 6 {m t , nib, it%h} in lattice units fulfill 

rh < 0.5 and rh ■ L Sit > 2, (46) 

which already is a rather loose condition in comparison with the common situation in QCD, where 
one usually demands at least rh ■ L s ,t > 3. In this model, however, the presence of massless Goldstone 
modes is known to induce algebraic finite size effects, which is why it is not meaningful to impose 
a much stronger constraint in Eq. (|46[) . since the quantity rh ■ L Stt only controls the strength of the 
exponentially suppressed finite size effects caused by the massive particles. 

Employing a top mass of 175 GeV and a Higgs boson mass of roughly 700 GeV, which will turn 
out to be justified after the upper mass bound has eventually been established, it should therefore be 
feasible to reach energy scales between 1400 GeV and 2800 GeV on a 32 4 -lattice. 

For the purpose of investigating the cutoff dependence of the upper mass bound a series of direct 
Monte-Carlo calculations has been performed with varying hopping parameters k associated to cutoffs 
covering approximately the given range of reachable energy scales. At each value of k the Monte-Carlo 
computation has been rerun on several lattice sizes to examine the respective strength of the finite 
volume effects, ultimately allowing for the infinite volume extrapolation of the obtained lattice results. 
In addition, a corresponding series of Monte-Carlo calculations has been performed in the pure <& 4 - 
theory, which will finally allow to address the question for the fermionic contributions to the upper 
Higgs boson mass bound. The model parameters underlying these two series of lattice calculations 
are presented in Tab. IIVI 
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TABLE IV: The model parameters of the Monte-Carlo runs underlying the subsequent lattice calculation of 
the upper Higgs boson mass bound are presented. In total, a number of 45 Monte-Carlo runs have been 
performed for that purpose. The available statistics of generated field configurations Nconf varies depending 
on the respective lattice volume. In detail we have Nconf ~ 20, 000 for 12 < L s < 16, Nconf ~ 10, 000— 15, 000 
for L s = 20, Nconf ~ 8, 000 - 16, 000 for L s = 24, N Co nf « 3, 000 - 5, 000 for L s = 32. The numerically 
determined values of 1/v and A are also approximately given. These numbers vary, of course, depending on 
the respective lattice volumes and serve here only for the purpose of a rough orientation. The degenerate 
Yukawa coupling constants in the upper four rows have been chosen according to the tree-level relation in 
Eq. (|35[) aiming at the reproduction of the phenomenologically known top quark mass. In the other cases it 
is exactly set to zero recovering the pure <E> 4 -tlieory. 



However, before discussing the obtained lattice results, it is worthwhile to recall what behaviour 
of the considered observables is to be expected from the knowledge of earlier lattice investigations. 
For the pure $ -theory and neglecting any double-logarithmic contributions the cutoff dependence 
of the Higg s boson mass as well as the renormalized quartic coupling constant has been found in 
Refs. @, El HI to be of the form 

= A m .[log(A 2 / /J 2 ) + B m ]- 1/2 , (47) 
A r = A\ ■ [log(A 2 //i 2 ) + B\\ 1 , (48) 

where fj, denotes some unspecified scale, and A m ^\ = A m> \{ji), B m \ = B mi \(/j,) are constants. Since 
this result has been established in the pure $ 4 -theory, it is thus worthwhile to ask whether these scaling 
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laws still hold in the considered Higgs- Yukawa model including the coupling to the fermions. In that 
respect it is remarked that the same functional dependence has also been observed in other analytical 
studies, for instance in Ref. [13] examining a Higgs- Yukawa model in continuous Euclidean space-time 
based, however, on an one-component Higgs field. In that study the running of the renormalizcd 
coupling constants with varying cutoff has been investigated by means of renormalized perturbation 
theory in the large A^-limit. Furthermore, the scaling behaviour of the renormalized Yukawa coupling 
constant has also been derived. It was found to be 

y r = A y -[log(A 2 /n 2 )+B y ]- 1/2 , (49) 

where A y = A y (fi) and B y = B y (fi) are again so far unspecified constants and y r stands here for 
the renormalized top and bottom Yukawa coupling constants yt, r and yb,r, respectively, as defined in 
Eq. 4M]). 

Now, the numerically obtained Higgs boson masses rriHp resulting in the above specified lattice 
calculations are finally presented in Fig. [BJ where panel (a) refers to the full Higgs- Yukawa model 
while panel (b) displays the corresponding results in the pure $ 4 -theory. To illustrate the influence 
of the finite lattice volume those results, belonging to the same parameter sets, differing only in the 
underlying lattice size, are connected by dotted lines to guide the eye. From these findings one learns 
that the model indeed exhibits strong finite volume effects when approaching the upper limit of the 
defined interval of reachable cutoffs, as expected. 
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FIG. 6: The Higgs propagator mass rriHp is presented in units of the vacuum expectation value v versus 1/v. 
These results have been determined in the direct Monte-Carlo calculations specified in Tab. IIVI Those runs 
with identical parameter sets differing only in the underlying lattice volume are connected via dotted lines to 
illustrate the effects of the finite volume. The dashed curves depict the fits of the lattice results according to 
the finite size expectation in Eq. (|52|) as explained in the main text. Panel (a) refers to the full Higgs- Yukawa 
model, while panel (b) shows the corresponding results of the pure "I> 4 -theory. 



In Fig. [6k, one sees that the Higgs boson mass seems to increase with the cutoff A on the smaller 
lattice sizes. This, however, is only a finite volume effect. On the larger lattices the Higgs boson mass 
decreases with growing A as expected from the triviality property of the Higgs sector. In comparison 
to the results obtained in the pure <J> 4 -theory shown in Fig. |B]d the aforementioned finite size effects, 
being of order 10 % here, are much stronger and can thus be ascribed to the influence of the coupling 
to the fermions. This effect directly arises from the top quark being the lightest physical particle in 
the here considered scenario. 

At this point it is worthwhile to ask whether the observed finite volume effects can also be under- 
stood by some quantitative consideration. For the weakly interacting regime this could be achieved, 
for instance, by computing the constraint effective potential 50, 51] (CEP) in terms of the bare model 
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parameters for a given finite volume as discussed in Ref. [38[ , which then allowed to predict the numer- 
ical lattice data for given bare model parameters. In contrast to that scenario the same calculation is 
not directly useful in the present situation, since the underlying (bare) perturbative expansion would 
break down due to the bare quartic coupling constant being infinite here. This problem can be cured 
by parametrizing the four-point interaction in terms of the renormalized quartic coupling constant. 
Starting from the definition of A r in Eq. (I31[) one can directly derive an estimate for the Higgs boson 
mass in terms of X r according to 
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which respects all contributions of order 0(\ r ). It is remarked that the above contribution Uf[v] 
contains all fcrmionic loops in the background of a constant scalar field and has already been discussed 
in Ref. [38] , while the underlying definition of the eigenvalues (p) of the free overlap operator with 
±Im(i/ ± (p)) > has been taken from Ref. (35| . 

Combining the above result with the expected scaling law given in Eq. (|4"5)l a crude estimate for 
the observed behaviour of the Higgs boson mass presented in Fig. [6] can be established according to 



<-He 



log(v 



v dv 



(52) 



where double- logarithmic terms have been neglected and A' x , B' x are so far unspecified parameters. 

Since the value of the renormalized quartic coupling constant A,, is not known apriori, the idea is 
here to use the result in Eq. ([52]) as a fit ansatz with the free fit parameters A' x and B' x to fit the 
observed finite volume behaviour of the Higgs boson mass presented in Fig. [5] These lattice data have 
been given in units of the vacuum expectation value v, plotted versus 1/v, to allow for the intended 
direct comparison with the analytical finite volume expression in Eq. (|52[) . The resulting fits are 
depicted by the dashed curves in Fig. |6j where the free parameters A' x and B' x have independently 
been adjusted for every presented series of constant lattice volume in the full Higgs- Yukawa model 
and the pure $ 4 -theory, respectively. Applying the above fit ansatz simultaneously to all available 
data does not lead to satisfactory results, since the renormalized quartic coupling constant itself also 
depends significantly on the underlying lattice volume, as will be seen later in this section. 

One can then observe in Fig. [6] that this fit approach can describe the actual finite volume cutoff 
dependence of the presented Higgs boson mass satisfactorily well, unless the vacuum expectation value 
v becomes too small. In that case the model does no longer exhibit the expected (infinite volume) 
critical behaviour in Eq. (I47H48[) which the derivation of the above fit ansatz was built upon. Staying 
away from that regime, however, the observed finite volume behaviour of the Higgs boson mass can 
be well understood by means of the analytical expression in Eq. (|52p . 

To obtain the desired upper Higgs boson mass bounds m^f(A) these finite volume results have to 
be extrapolated to the infinite volume limit and the renormalization factor Zq has to be properly 
considered. For that purpose the finite volume dependence of the Monte-Carlo results on the renor- 
malized vev v r = v I \J Zq and the Higgs boson mass rriHp, as obtained for the two scenarios of the 
full Higgs- Yukawa model and the pure $ 4 -theory, is explicitly shown in Fig. [7] One sees in these 
plots that the finite volume effects are rather mild at the largest investigated hopping parameters k 
corresponding to the lowest considered values of the cutoff A, while the renormalized vev as well as the 
Higgs boson mass itself vary strongly with increasing lattice size L s at the smaller presented hopping 
parameters, as expected. 

It is well known from lattice investigations of the pure $ 4 -theory [44Tj46| that the vev as well as 
the mass receive strong contributions from the Goldstone modes, inducing finite volume effects of 
algebraic form starting at order 0(Lj 2 ). The next non-trivial finite volume contribution was shown 
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to be of order 0(L S 4 ). In Fig. [7] the obtained data are therefore plotted versus l/L 2 . Moreover, the 
aforementioned observation justifies to apply the linear fit ansatz 

fV) n (L- 2 ) = AW m +B2] m -L- 2 (53) 

to extrapolate these data to the infinite volume limit, where the free fitting parameters Av} m an d 
Bv}m with the subscripts v and m refer to the renormalizcd vcv v r and the Higgs boson mass rriH p , 
respectively. 
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FIG. 7: The dependence of the renormalized vev v r = v/s/Zg and the Higgs propagator mass map on the 
squared inverse lattice side length 1/L 2 g is presented in the upper and the lower panel rows, respectively, as 
determined in the direct Monte-Carlo calculations specified in Tab. IIVI Panels (a) and (c) show the results 
for the full Higgs- Yukawa model, while panels (b) and (d) refer to the pure <J? 4 -theory. In all plots the dashed 
curves display the parabolic fits according to the fit ansatz in Eq. (|54[) , while the solid lines depict the linear 
fits resulting from Eq. (|53p for the two lower threshold values L' 3 = 16 (red) and L' s — 20 (black). 



To take the presence of higher order terms in 1 / L 2 S into account only the largest lattice sizes are 
included into this linear fit. Here, we select all lattice volumes with L s > L' s . As a consistency 
check, testing the dependence of the resulting infinite volume extrapolations on the choice of the 
fit procedure, the lower threshold value L' s is varied. The respective results are listed in Tab. [V] 
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TABLE V: The results of the infinite volume extrapolations of the Monte-Carlo data of the renormalized vev 
v r and the Higgs boson mass rriH-p are presented as obtained from the parabolic ansatz in Eq. ([54} and the 
linear approach in Eq. (|53[) for the considered lower threshold values L' s — 16 and L' s — 20. The final results 
on v r and mHp, displayed in the very right column, are determined here by averaging over the parabolic and 
the two linear fit approaches. An additional, systematic uncertainty of these final results is specified in the 
second pair of brackets taken from the largest observed deviation among all respective fit results. 

Moreover, the parabolic fit ansatz 

f(p) (L~ 2 ) = A {p) + B {p) -L- 2 + C {p) ■L~ i (54) 

Jv,m\ s I v,m 1 v,m s 1 ^v.rn s V" / 

is additionally considered. It is applied to the whole range of available lattice sizes. The deviations 
between the various fitting procedures with respect to the resulting infinite volume extrapolations 
of the considered observables can then be considered as an additional, systematic uncertainty of the 
obtained values. The respective fit curves are displayed in Fig.[7]and the corresponding infinite volume 
extrapolations of the renormalized vev and the Higgs boson mass, which have been obtained as the 
average over all presented fit results, are listed in Tab. [V] 

The sought-after cutoff-dependent upper Higgs boson mass bound, and thus the main result of 
this paper, already presented in Fig. [TJ can then directly be obtained from the latter infinite volume 
extrapolation. The bounds arising in the full Higgs- Yukawa model and the pure <I> 4 -theory are jointly 
presented in Fig. [T^,. In both cases one clearly observes the expected decrease of the upper Higgs 
boson mass bound with rising cutoff A. Moreover, the obtained results can very well be fitted with 
the expected cutoff dependence given in Eq. (|4"T)) . as depicted by the dashed and solid curves in Fig.QJi, 
where A m , B m are the respective free fit parameters. 

Concerning the effect of the fermion dynamics on the upper Higgs boson mass bound one finds 
in Fig. QJ, that the individual results on the Higgs boson mass in the full Higgs- Yukawa model and 
the pure $ 4 -theory at single cutoff values A are not clearly distinguishable from each other with 
respect to the associated uncertainties. Respecting all presented data simultaneously by considering 
the aforementioned fit curves also does not lead to a much clearer picture, as can be observed in 
Fig. [1^, where the uncertainties of the respective fit curves are indicated by the highlighted bands. At 
most, one can infer a mild indication from the presented results, being that the inclusion of the fermion 
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dynamics causes a somewhat steeper descent of the upper Higgs boson mass bound with increasing 
cutoff A. A definite answer regarding the latter effect, however, remains missing here due to the size 
of the statistical uncertainties. The clarification of this issue would require the consideration of higher 
statistics as well as the evaluation of more lattice volumes to improve the reliability of the above 
infinite volume extrapolations. 

On the basis of the latter fit results one can extrapolate the presented fit curves to very large values 
of the cutoff A as illustrated in Fig. [TJd. It is intriguing to compare these large cutoff extrapolations 
to the results arising from the perturbative consideration of the Landau pole presented, for instance, 
in Ref. 19]. One observes good agreement with that perturbatively obtained upper mass bound even 
though the data presented here have been calculated in the mass degenerate case and for Nj = 1. 
This, however, is not too surprising according to the observed relatively mild dependence of the upper 
mass bound on the fermion dynamics. 

For clarification it is remarked that a direct quantitative comparison between the aforementioned 
perturbative and numerical results has been avoided here due to the different underlying regularization 
schemes. With growing values of A, however, the cutoff dependence becomes less prominent, thus 
rendering such a direct comparison increasingly reasonable in that limit. 




(a) (b) (c) 

FIG. 8: The dependence of the renormalized quartic coupling constant A r as well as the top quark mass mt on 
the squared inverse lattice side length 1/L 2 S is presented as calculated in the direct Monte-Carlo calculations 
specified in Tab. IIVI Panels (a) and (c) show the results for the full Higgs- Yukawa model, while panel (b) 
refers to the pure $ 4 -theory. In all plots the dashed curves display the parabolic fits according to the fit ansatz 
in Eq. (|54p . while the solid lines depict the linear fits resulting from Eq. (|53[) for the two lower threshold values 
L' s = 16 (red) and L' a = 20 (black). 

Furthermore, the question for the cutoff dependence of the renormalized quartic coupling constant 
A r and - in the case of the full Higgs- Yukawa model - the top quark mass with its associated value 
of the renormalized Yukawa coupling constant yt lT shall be addressed. For that purpose we follow 
exactly the same steps as above. The underlying finite volume lattice results on the renormalized 
quartic coupling constant and the top quark mass are fitted again with the parabolic and the linear 
fit approaches in Eq. (|53|) and Eq. (|54|) as presented in Fig. [8j 

The corresponding infinite volume extrapolations are listed in Tab. I VI [ where the final extrapolation 
result is obtained by averaging over all performed fit approaches. An additional systematic error is 
again estimated from the deviations between the various fit procedures. 

The sought-after cutoff dependence of the aforementioned renormalized coupling constants can then 
directly be obtained from the latter infinite volume extrapolations. The respective results are presented 
in Fig. [5] and within the achieved accuracy one observes the renormalized coupling parameters to 
be consistent with the expected decline when increasing the cutoff A as expected in a trivial theory. 
Again, the obtained numerical results arc fitted with the analytically expected scaling behaviour given 
in Eq. (|48[) and Eq. (fi^)) . As already discussed for the case of the Higgs boson mass determination, 
the individual measurements of X r in the two considered models at single cutoff values A are not 
clearly distinguishable. Respecting the available data simultaneously by means of the aforementioned 
fit procedures also leads at most to the mild indication that the inclusion of the fermion dynamics 
results in a somewhat steeper descent of the renormalized quartic coupling constant with rising cutoff 
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TABLE VI: The results of the infinite volume extrapolations of the Monte-Carlo data of the renormalized 
quartic coupling constant A r and the top quark mass mt are presented as obtained from the parabolic ansatz 
in Eq. (154[) and the linear approach in Eq. (|53f) for the considered lower threshold values L' B — 16 and L' B — 20. 
The final results on A r and mt, displayed in the very right column, are determined here by averaging over 
the parabolic and the two linear fit approaches. An additional, systematic uncertainty of these final results 
is specified in the second pair of brackets taken from the largest observed deviation among all respective fit 
results. 

A as compared to the pure $ 4 -theory. A definite conclusion in this matter, however, cannot be drawn 
at this point due to the statistical uncertainties encountered in Fig. |H1 

Finally, the renormalized Yukawa coupling constant is compared to its bare counterpart depicted 
by the horizontal line in Fig. O). Since the latter bare quantity was chosen according to the tree-level 
relation in Eq. (1551) aiming at the reproduction of the physical top quark mass, one can directly infer 
from this presentation how much the actually measured top quark mass differs from its targeted value 
of 175 GeV. Here, one observes a significant discrepancy of up to 2 %, which can in principle be fixed 
in follow-up lattice calculations, if desired. According to the observed rather weak dependence of 
the upper Higgs boson mass bound on the Yukawa coupling constants, however, such an adjustment 
would not even be resolvable with the here achieved accuracy. 



VIII. SUMMARY AND CONCLUSIONS 

The aim of the present work has been the non-perturbative determination of the cutoff-dependent 
upper mass bound of the Standard Model Higgs boson based on first principle computations, in 
particular not relying on additional information such as the triviality property of the Higgs- Yukawa 
sector. The motivation for the consideration of the aforementioned mass bound finally lies in the 
ability of drawing conclusions on the energy scale A at which a new, so far unspecified theory of 
elementary particles definitely has to substitute the Standard Model, once the Higgs boson and its 
mass m,H will have been discovered experimentally. In that case the latter scale A can be deduced by 
requiring consistency between the observed mass mu and the upper and lower mass bounds ui^f(A) 
and m 1 ™ (A) , intrinsically arising from the Standard Model under the assumption of being valid up 
to the cutoff scale A. 

The Higgs boson might, however, very well not exist at all, especially since the Higgs sector can only 
be considered as an effective theory of some so far undiscovered, extended theory, due to its triviality 
property. In such a scenario, a conclusion about the validity of the Standard Model can nevertheless 
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FIG. 9: The cutoff dependence of the renormalized quartic and Yukawa coupling constants is presented in 
panels (a) and (b), respectively, as obtained from the infinite volume extrapolation results in Tab. IVll The 
dashed and solid curves are fits with the respective analytically expected cutoff dependence in Eq. (|48l) and 
Eq. (|49|) . The horizontal line in panel (b) indicates the bare degenerate Yukawa coupling constant underlying 
the performed lattice calculations. 

be drawn, since the non-observation of the Higgs boson at the LHC would eventually exclude its 
existence at energies below, lets say, 1 TeV thanks to the large accessible energy scales at the LHC. 
An even heavier Higgs boson is, however, definitely excluded without the Standard Model becoming 
inconsistent with itself according to the results in section I VIII and the requirement that the cutoff A 
be clearly larger than the mass spectrum described by that theory. In the case of non-observing the 
Higgs boson at the LHC after having explored its whole energy range, one can thus conclude on the 
basis of the latter results, that new physics must set in already at the TeV-scale. 

For the purpose of establishing the aforementioned cutoff-dependent mass bound, the lattice ap- 
proach has been employed to allow for a non-perturbative investigation of a Higgs- Yukawa model 
serving as a reasonable simplification of the full Standard Model, containing only those fields and 
interactions which are most essential for the Higgs boson mass determination. This model has been 
constructed on the basis of Liischer's proposals in Ref. [3(| for the construction of chirally invari- 
ant lattice Higgs- Yukawa models adapted, however, to the situation of the actual Standard Model 
Higgs-fermion coupling structure, i.e. for tp being a complex doublet equivalent to one Higgs and 
three Goldstone modes. The resulting chirally invariant lattice Higgs- Yukawa model, constructed 
here on the basis of the Neuberger overlap operator, then obeys a global S , [/(2)l x U(1)y symmetry, 
as desired. 

The fundamental strategy underlying the determination of the cutoff-dependent upper Higgs boson 
mass bounds has then been the numerical evaluation of the maximal Higgs boson mass attainable 
within the considered Higgs- Yukawa model in consistency with phenomenology. The latter condition 
refers here to the requirement of reproducing the phenomenologically known values of the top and 
bottom quark masses as well as the renormalized vacuum expectation value v r of the scalar field, 
where the latter condition was used here to fix the physical scale of the performed lattice calculations. 
Owing to the potential existence of a fluctuating complex phase in the non-degenerate case, the top 
and bottom quark masses have, however, been assumed to be degenerate in this work. Applying this 
strategy requires the evaluation of the model to be performed in the broken phase, but close to a 
second order phase transition to a symmetric phase, in order to allow for the adjustment of arbitrarily 
large cutoff scales, at least from a conceptual point of view. 

As a first step it has explicitly been confirmed by direct lattice calculations that the largest attain- 
able Higgs boson masses are indeed observed in the case of an infinite bare quartic coupling constant, 
as suggested by perturbation theory. Consequently, the search for the upper Higgs boson mass bound 
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has subsequently been constrained to the bare parameter setting A = oo. The resulting finite vol- 
ume lattice data on the Higgs boson mass turned out to be sufficiently precise to allow for their 
reliable infinite volume extrapolation, yielding then a cutoff-dependent upper bound of approximately 
niff(A) — 630 GeV at a cutoff of A = 1500 GeV. These results were moreover precise enough to actu- 
ally resolve their cutoff dependence as demonstrated in Fig. [TJ which is in very good agreement with 
the analytically expected logarithmic decline, and thus with the triviality picture of the Higgs- Yukawa 
sector. 

It is remarked here, that this achievement has been numerically demanding, since the latter loga- 
rithmic decline of the upper bound rn^f (A) is actually only induced by subleading logarithmic contri- 
butions to the scaling behaviour of the considered model close to its phase transition, which had to be 
resolved with sufficient accuracy. By virtue of the analytically expected functional form of the cutoff- 
dependent upper mass bound, which was used to fit the obtained numerical data, an extrapolation 
of the latter results to much higher energy scales could also be established, being in good agreement 
with the corresponding perturbatively obtained bounds [19]. A direct comparison has, however, been 
avoided due to the different underlying regularization schemes. 

The interesting question for the fermionic contribution to the observed upper Higgs boson mass 
bound has then been addressed by explicitly comparing the latter findings to the corresponding results 
arising in the pure $ 4 -theory. For the considered energy scales this potential effect, however, turned 
out to be not very well resolvable with the available accuracy of the lattice data. The performed 
fits with the expected analytical form of the cutoff dependence only mildly indicate the upper mass 
bound in the full Higgs- Yukawa model to decline somewhat steeper with growing cutoffs than the 
corresponding results in the pure $ 4 - theory. To obtain a clearer picture in this respect, higher accuracy 
of the numerical data and thus higher statistics of the underlying field configurations would be needed. 
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